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■^ ' Abstract 

In previous work [12] the authors found integral formulas for probabihties 

in the asymmetric simple exclusion process (ASEP) on the integer lattice Z. 
fT^ . The dynamics are uniquely determined once the initial state is specified. In this 

^ I note we restrict our attention to the case of step initial condition with particles 

at the positive integers Z+ and consider the distribution function for the mth 
^f) ■ particle from the left. In [12] an infinite series of multiple integrals was derived 

for the distribution. In this note we show that the series can be summed to 
^f ' give a single integral whose integrand involves a Fredholm determinant. We 

use this determinant representation to derive (non-rigorously, at this writing) 






o 

OO 



f^ , a scaling limit. 



% 



I. Introduction 



The asymmetric simple exclusion process (ASEP) is a basic interacting particle 
model for nonequilibrium phenomena. Since its introduction by Spitzer [10] in 1970, 
it has become a popular and much studied model. (See, e.g. [3llll[6],|9].) In ASEP on 
the integer lattice Z particles move according to two rules: (1) A particle at x waits 
an exponential time with parameter one (independently of all other particles), and 
then it chooses y with probability p{x,y); (2) If y is vacant at that time it moves to 
y, while if y is occupied it remains at x and restarts its clock. The adjective "simple" 
refers to the fact that allowed jumps are one step to the right, p{x,x + 1) = p, or 
one step to the left, p{x, x — 1) = 1— p = q. Because we deal with continuous time, 
we need not worry about two or more particles attempting jumps at the same time. 



The model is called T(totally)ASEP if either p = 1 (particles hop only to the right) 
or g = 1 (particles hop only to the left). 

The dynamics are uniquely determined once we specify the initial state. We 
restrict our attention to the case of step initial condition with particles at the positive 
integers Z"*". With this initial condition it makes sense to talk about the position of 
the mth left-most particle at time t. We denote this (random) position by Xmit). (So 
Xm{0) = m.) 

In a now classic paper, Johansson [2] for the case of TASEP with step initial 
condition showed that the distribution of Xm(t) is related to the distribution of the 
largest eigenvalue in the unitary Laguerre ensemble of random matrix theory. This 
connection with random matrix theory leads to an expression for the distribution 
of Xm{t) as a Fredholm determinant. The importance of the Fredholm determinant 
representation is that it makes possible an analysis of the regime of universal fluctua- 
tions. (See Johansson's Corollary 1.7. For further discussion of this universal regime 
see pin].) 

For the general ASEP model on Z there is, as far as the authors know, no known 
connection to random matrix theory and no analysis of the analogous scaling regime 
as analyzed by Johansson in the TASEP case. (However for stationary ASEP Balazs 
and Seppalainen jT| and Quastel and Valko [7] prove that the variance of the current 
across a characteristic is of order t^/'^ and that the diffusivity has order t^/^.) In recent 
work [12] we showed (building on ideas from Bethe Ansatz [8]) that F{xm(t) < x) 
can be expressed as an infinite series where the kth order term is a /c-dimensional 
integral. 

In this paper we show that when p and q are nonzero this infinite series can be 
summed to give a single integral whose integrand involves a Fredholm determinant. 
We use this determinant representation to derive a scaling limit for fixed m with 
X, t —>■ oo. (It is not an extension of the scaling limit of Johansson since m is fixed.) 
It remains conjectural since the derivation lacks a justification for an interchange of 
limits, but it is surely true. 

II. Determinant Representation 

The result is stated in terms of an operator with kernel 

p + q^^ - C 

where 

If g 7^ and Cr denotes the circle with center zero and radius R then K{C,, C,') is a 
smooth function on Cr x Cr when R is sufficiently large. Then the operator K on 



L\Cr) defined bjfl 

Kf{0= f K{^,0f{0de 
JCr 

is trace class. With the notations 

r = p/g, (A; r)^ = (1 - A) (1 - Ar) ■ ■ ■ (1 - A r^'^), 
the result is that when p and q are both nonzero, 

dei{I - XqK) d\ 



^{Xm{t)<x) 



(A;r) 



A 



where the integral is taken over a circle so large that all the singularities of the 
integrand lie inside it. Evaluating the integral by residues gives the equivalent formula 



m— 1 






det {I -qr-' K) 



-rJ- 



J < m ~ 1 



(When r = 1 this is modified in the obvious way.) In particular, 

P(xi(t) > x) = det{I -qK). 



To state the Corollary to Theorem 5.2 of [12], from which ([T]) will follow, we recall 
the the definition of the r-binomial coefficient 



n 



l-T^)(l-r 



A^-l^ 



'1-r 



Af-n+l^ 



(l-r)(l-r2)...(l-r") 
The result was that if the initial state is Z+ and q ^ theru 



{Xm{t) =X) = (-1)™+^ Y. h 



k>m 



k-1 

k — m 



{k~m){k-m+l)/2 k'm+{k~m){k+m-l)/2 



n 



- 6 



1-6 ■■■a 



Icn Jcn 14 P + ?66 - 6 J](l - e,) (g6 - P) 

i 

where in the fc-dimensional integral all indices run over {1, . . . , k}. 



n(erv(«'^*)rfei---rfa, 



^AU contour integrals are to be given a factor l/27ri. 

^The coefficients were given in |12| in terms of what we denoted there by 



N 



which is related 



to 



N 



by 



N 



^ qUiN-n) 



N 



If we sum this over x from — oo to a; (which we may do since we may take R> 1] 
we obtain 



[X, 



.{t) <x) = i-ir E ^ 



X 



Cr 



k>m 
- ^^ 



Cn t/j P + ^^^^i " ^* 



k-1 

k — m 



(fc-m)(fc-m+l)/2 km+{k-m)(k+m-l)/2 
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- 6) (Q^i - P) \ 



nfe^e^^^'^V^i-'-^e^- (2) 



The observation that will allow us to express this in terms of Fredholm determi- 
nants is the identity 



det 



1 



p + Qiiij - ii J I 



l<i,j<k 



.l)k^pq)Hk-i)/2Y^ 



e, - e. 



f_i^ p + qiiij - ii 



i.Ml-£,)fQ 



- it){.(iii - p) 

(3) 
This can be seen as follows. We may assume t 7^ 1 (i.e., q 7^ 1/2) since both sides 
are continuous in q. If we make the substitutions 



6 



r?i + l 



-1 ' 



then 



Since 



1 



Tli + T" 

1 (l + rr7i)(l + rr7,) 



p + qiiij-ii p(i-^) 

det ' 



Vi - ^Vj 



(4) 



is a Cauchy determinant we can evaluate it, and we find that the determinant of (HI) 
equals 



*7^i 



^i - ^^j 



We compute 



Vi 
Vi - Vj 



p(l - r) 



? 



(ge,-p)(l-e, 
e. - e. 



r/i - rr^j- p + q^i^j - ^i 

and find that ([5]) equals the right side of (IHl). 

From ([3]) we see that when also p 7^ ([2]) may be written 



p(x^(t) <x) = i-ir Yl 



k>in 



k\ 



k 
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rn(m— 1)/2— fc(m— 1) — m(»n— l)/2+fcm 



X / • • ■ / det{K{^i, Q)i<i,j<k d^i--- dC,k 
'Cr Jcr 



C_x)™ T-'^(™-i)/2 y^ 
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The last integral is a coefficient in the Fredholm expansion of det(/ — \K). In 



fact 



E 

fc=0 



(-A)* 



k\ 



Cr J Cr 



det(ir(e„ ij))x<i,J<^^ di^--- d^k = det(/ - XK). (7) 



From [5^, P-26] we have for \z\ small enough 



k>m 



k-l' 

k — m 



z^ = z'^ 



E 

i>o 



m + j — 1 
J 



z^ = z"' 



m— 1 ^ m 

n^^=n — ^ 

J-J- l-T'Z -'--'- 1-T'^-^Z 
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If we set z = pr ™ A this gives for | A | large enough 



fl\m ^m{m~l)/2 V^ 



k>m '- 



k-l 

k — m 



m ^ 






If we multiply ([7]) and ([H]), multiply by A ^, and integrate and we see that (^ is the 
same as ([1]). 

Remark. Formula ([1]) holds when p and g are nonzero. When p ^ we must get 
the TASEP determinant for the probability, and the question arises whether this is 
easy to see. The answer seems to be that it is not easy, but it can be derived from it 
with some work. 

III. A Scaling Conjecture 

Denote by Kq the operator on L^(]R) with kernel 



Ko{z, z') 



(p^W) {z'^+z'^)/A+pqzz' 



^27r 



(This is the symmetrization of the Mehler kernel.) When p ^ q this is trace class. 
The conjecture is that if p < g, so there is a drift to the left, we have for each m 

hm P (x^(t) < (p - q)t + (g - p)yt^l-) = [ ^''^' " f ^° ^^'^-^^ ^, (9) 

where, as before, the integral is taken over a circle so large that all the singularities 
of the integrand lie inside it. 
We know from ([T]) that 

<lei{I - XqK) d\ 



F (x„(t) <{p- q)t + {q- p)yt'/^) = f 



(A; T)m A ' 
where K is as before and 

x= (p-q)t + {q-p)yf/'^. (10) 

Therefore the conjecture would follow if we can show that 
det(/ — XK) -^ det(/ — A-K'o^(-j/,oo)) as t — > oo, uniformly on compact A-sets. (11) 

The Fredholm determinants are entire functions of A, and the coefficients in their 
expansions about A = are universal polynomials in the traces of powers of the 
operators. By considering the coefficients successively we can see that if flTTl) is true 
then necessarily 

(i) trK" ^tT{KoX(^^y^oo))T ast^oo for n= 1,2, ■■■. 

If (i) holds then each coefficient in the expansion of det(/ — XK) converges to 
the corresponding coefficient in the expansion of det(/ — A/<'oX(_j;,oo))- This is not 
sufficient to give fITT]) . It would be sufficient if we also knew that 

(ii) det(/ — XK) is uniformly bounded for large t on compact X-sets. 

For if this holds one sees from the Cauchy inequalities that the convergence of the 
series for det(/ — XK) is uniform in t. 

We shall show that (i) holds. 

First, though, we remark that if the conjecture is true then we expect the right 
side of dn]) to be a distribution function in y for each m. As ?/ — *■ — cxo the numerator 
in the integrand approaches one, and expanding the contour shows that the resulting 
integral equals zero. But what about the y -^ +oo limit, which should equal one? 
The right side becomes 

^ "^ detjl-qr-' Kq) 

''h n (1 - r'-) • <''' 

J < m — 1 



where Kq acts on all of M. If this is to equal one when m = 1 then we must have 
det(/ — qKo) = 0. If this holds and flT2|) is equal to one when m = 2 then we must 
also have det(J — qT'^Ko) = 0. And so on. The conclusion is that det (/ — gr~* Kq) 
should equal zero for all i > 0, in other words that all r^q should be eigenvalues of 
Kq. These are indeed eigenvalues and the corresponding eigenfunctions are 



^-i,-~p-)zV4^ I jtjZt. 



Z 



where the Hi are the Hermite polynomials. 

Next, we explain where the conjecture came from: there is a kernel with the same 
Fredholm determinanio as K that converges pointwise to a kernel with the same 
Fredholm determinant as -K'o^(-3/,oo)- 

When p < q we have, when |^| and |^'| are large enough. 

Therefore if we define 

A(e, z) = {q-pf- e-^^^, B{z, O = — e'^^^?-^, 

qi-p qi' -p 



then 



mA')= / A{i,z)B{z,Odz. 
Jo 



The kernel 

K,{z, z') = [ Biz, OA(e, z')d^ = ^ [ f-^e^^^-^-'^W^d^ (13) 
Jcn 27rz Jcj, [q^ - pf 

on R+ has the same Fredholm determinant as K on Cr. 
The part of the exponent with the factor t, 

(p-g)loge + £(0, (14) 

has a critical point at ^ = 1, its second derivative is positive there, and on the line 
Re^ = 1 its real part has an absolute maximum there. When x < we may replace 
Cr by this line, and we make the substitution ^ — * 1 + i^t~^/^. We also make the 
substitutions 

z^{q- p)zt^/^, z' ^ (q- p)z't^'^, 



•^We use the term "Fredholm determinant" to mean the infinite series. In case the kernel is 
continuous and trace class on some Hilbcrt space this is the same as the operator determinant. 
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and multiply the kernel by (g — p) t^/^, so the Fredholm determinant is unchanged. 
We let t -^ oo and obtain the pointwise limit 

K2{z, z') = ^e-('?"-P"'-(^-P)^)'/2_ 



Here z, z' G M^. We may replace this by the kernel -^'3(2;, z')X(-y^oo){z')i where 
z, z' eM. and 

K^{z, z') = J—e-^<i--v''fl'^, 

The kernel Kq{z^ z') is the symmetrization of K^{z, z'). In fact 
Ko{z, z') = e^'^-P'^^'/^Ksiz, y) e-(«-P)^"/^ 

so they have the same Fredholm determinants. Hence, the conjecture. 

Finally, we establish (i). Instead of K we may use the kernel Ki given by flT3l) 
because the traces of the powers are the same. We make the variable change 

V 



and find that with x given by (11 01) 



27r2 J 



where 

f[n) = {q- p)y log ( 

\l-qr] 

g(n) = (p-q) log I ) + (g - pf 



1 — qrj J {1 — prj) {1 — qrj)' 

and 7 is a little circle around 1/q described clockwise. (The apparent complication 
here will make things simpler later.) 

The function g{ri) equals flT^ after our substitution. The line Re^ = 1 corresponds 
to the circle F with diameter (0, 1/q)- On this contour g has a critical point at r/ = 0, 
its second derivative is positive there, and the real part of g{r]) has an absolute 
maximum there. 



We consider first 



tTK, = ^—— /'e^('')*'''+^(^)*^. (15) 

q-p2-ni J^ 7] 



In the neighborhood of rj = 

fiv) = {q- Pfyn + 0{rf), gir]) = {q- pf rf /2 + 0{r^''). 

Since the singularity at 1/p > 1/q is outside 7 and F we may deform 7 to F if we 
replace the denominator i] in (TT5|) by r]+0. We then make the substitution i] — ^ irft'^/"^ 
and find (recalling that 7, and so F, is described clockwise) that 

t^cc q-p27rj^^ 17] + 

„ „ fOO POO T /-oo 

2vr J_^ Jq v^ io 



V27r j-j; 
Now comes the tricky part. We compute that 

where we set r/n+i = Vi- Recalling that r = p/q we may write this as 

1 f f eJ:,{t'^'fM+t9M) 



We want to deform all contours 7 to F. Suppose we deform the r^j-contour 7, 
while perhaps some are still 7 and others are already F. Since r < 1 the pole at 
7]j = TTjjj^i will be crossed in the deformation whether r/j+i is on 7 or F, while the 
pole at T^^rjj^i will not be crossed. The residue at the crossed pole will give rise to 
a lower-order integral. In the deformation of all contours to F we will eventually get 
a sum of integrals over F of order < n. But there is the problem that the integrands 
in these integrals will be singular at zero, so we cannot blithely deform all contours 
to F, even if we keep track of the residues. 

What we do is first expand F slightly to a circle F' with diameter (—a, l/q)-, where 
a is very small and positive. Because F is expanded, but just a little, we will still 
cross poles at the Trjj (for tjj on the contour 7 or the entire contour F') but not at 
the T~^rjj. (If we had shrunk F instead we would cross the pole at T'qj only for some 
rjj G F', which would complicate determining the lower-order integrals.) 

Expanding all contours as described we get a sum of integrals over F' of order 
< n, one for each subset S* of {1, . . . ,n}, in which the integrands are nonsingular. 

9 



The integral corresponding to 5* is obtained as follows: for each j & S remove the 
factor rjj — ttjj^i from the denominator, make the replacement rjj — > Tr^^+i in the 
remainder of the integrand, and multiply by 27ri. The integral is taken with respect 
to the rjj with j ^ S. This is the integrand we get after we take the residues that 
arise from crossing the poles at the rjj = Trjj^i with j G Su 

The result may be described equivelently as follows: replace the product 

j 
with all rij G 7 by 



j 
with all rjj G F'. It follows that flT6l) is equal to 

Once we have this we can take the t ^ 00 limit. Recall that in the above a could 
have been any sufficiently small positive number. We take it to be t~^/^, so the left- 
most point of r' is — 1~^/^, and V is vertical there. On V we make the substitutions 
rjj — > irij t^^l"^ and obtain in the t ^ 00 limit 



-1 poo+i poo+i 



n 



J 



Q-{g-p?Vj /^+i{q-pfyr]3 



^ n V'^^ ~ ^'^i+i) ^ + 27rz(5(r7j - rr/j+i) rfr^i ■ ■ ■ dr]^ (17) 

j 
We used here the fact that 6 is homogeneous of degree —1. 

We now undo what we did before. Suppose we had the integral 

/00—i poo—i 

-00— i J— 00— i 



{2TTiqy 



and wanted to integrate over Imr^j = +1 instead over Imr^j = —1. If we raised 
the ?7j-contours successively we would obtain a sum of integrals of lower order over 
Imrij = 1 that come from the residues at various poles at rjj = Trjj^i. The sum of all 
integrals that arise is precisely equal to ( TT71) . Thus the limit of trK" is equal to (IT8|) . 



^Although 7 is expanded outward we multiply by 27rj rather than — 27ri because 7 is described 
clockwise. 

10 



We now substitute into flTSj) the integral representations, valid when Imrij = — 1, 

1 



e 



<V3~rr)j+i)zj ^^ 



and integrate first with with respect to the rjj. We obtain 



'J' 



-I /"OO /"OO "' 

i2nr/'Jo Jo U 



where zq = Zn- This is equal to 

1 /'OO /'OO ^ 

— — / . . . / n e-(^^— ^^^)^/^ dzi ■ ■ ■ rfz„ = tr (7^3 X(_,,oo))" = tr (i^o X(_,,oo))". 

l^TTJ J_j^ J_j^ _^.^_^ 

Thus (i) is established. What would it take to establish (ii)? It would be enough to 
find kernels that have the same Fredholm determinants as K and that have bounded 
trace norms on some Hilbert space. We could use the Hilbert-Schmidt norm instead 
(which is weaker and easier to compute) since that would give uniform boundedness 
of the regularized 2-determinant, and since the traces are bounded the determinant 
would also be. But so far we have not found any such kernels. 
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